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fON_ Objective 


Model propellant slosh for Europa Clipper 
using two pendulums such that controls 
engineers can predict slosh behavior 
during the mission. 


BACKGROUND 


* Importance of predicting propellant slosh 


— Sloshing changes CM (center of mass) of spacecraft and exerts 
forces and torques on spacecraft 


— Avoid natural frequencies of structures 
— Size ACS (Attitude Control Systems) thrusters to counteract forces 
and torques 
* Can model sloshing fluid as two pendulums with specific 
parameters (mass, length, damping) 


f, " Background 


+ Europa Clipper tanks and PMD 
— Bipropellant system 
— Cylindrical with domed top and bottom 
— 8-vane PMD (propellant management device) 
» CFD (computational fluid dynamics) data 
used as “real” slosh behavior 
— Have data for two propellants at three fill CED Simulation 
fractions each Solution Time 2 (s) 
— Initial condition of 15 degree free surface 
offset, released and allowed to settle 
— CFD requires long computing time -> Need a 
computationally simple model 


i, i Background 


+ Pendulum model 


— Model fluid movement as two pendulums Forces exerted on 
attached to central axis of the tank tank by fluid 
— For each CFD data set, find parameters: 
mass, frequency, damping ratio, attachment —mib\ 
height 
mL 


CM(t) = mLsiné@ (t) 


: fo. 
= mLsin Oe7St ( sin (ot —e t) + cos (ot —@ d) 
wf1-f 


d. \ Existing Literature 


+ SP-106 (1966), SwRI (2000): Sector tank mode Annular tank 
Analytical equations and empirical MoE ten; Smene (ein view) 
correlations for damping and PM 
frequency 

— Includes bare cylindrical (no PMD), Tank. 
sector, and annular tanks Wall 


Cassini paper illustration of 
+ Cassini slosh paper (1994): Two ‘double pendulum rmpite) 
pendulum model 
— Slosh around PMD was modeled as 
combination of sector and annular 
slosh modes 
— Two separate pendulums to model 
two slosh modes 
— Static mass component at bottom that 
experiences little movement 


METHODS OVERVIEW 


fan Generate CFD Data 


Solution Time 2 (s) 


* Propellants: NTO and MMH 
* Fill fractions: 25%, 50%, 85% 
+ Data: CM, Force, Moment (all 3 axes) 


Center of mass in xvs. time 


a 9 m0 10 1m 40 180 
Time () 


180200 


Curve fitting by finding 
parameters in 
pendulum equation that 
most closely match 
CFD 


Trying to resolve CFD 
into two pendulums 
Peak-to-peak values -> 
Initial guesses for 
damping and frequency 
of each pendulum 

Note much higher 
damping before first 
peak 


d. " Find Parameters to Fit CM Data 


Distance rm central axis 


Center of mass in x coordinate 


Matlab’s fsolve(x) -> 
Mass, damping, and 
frequency parameters to fit 
CMx CFD data 

Refine and iterate 


Distance rom canal acs 


Genter of mass in x coordinate 


lea Compare Sum of Pendulums to CFD Data 


10% (OM in Xv Time 


+ Sum of two pendulums 
generates model for 
propellant slosh 

* Should match both CM 
and Force data 


canny 


d. " Mean Error in Force Nr 


Force in X ve Time 


[= eraum | 


~* © 0 2 mM & 2 9 w 
Timees) 


+ Metric to quantify accuracy of fit: mean absolute difference 
between CFD force and pendulum model force 


Ve 
= > abs(CFD — pendulum) 


T 
* Select methods that minimize this 


RESULTS AND LITERATURE 
COMPARISON 


d. " Basis for results 


* Coordinate system — origin at top 
of tank 
+ Parameters prioritized fitting the 
behavior after the first peak 
+ Two pendulum model is an 
approximation only 
— PMD does not create a perfectly 
sector nor annular tank and is only a 
fraction of tank height 
— Parameters not constant over time 
— Model does not scale well with high 
fluid displacements 


Mass Participation Fraction NASA 


Mass participation fraction vs. fill fraction 
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Fill fraction 


+ Pendulum mass as a fraction of total fluid mass 
+ Monotonic trends 
+ Mass fractions are identical between NTO and MMH 
+ Piecewise linear fit 
— First two fill fractions — fluid partially submerges PMD, sloshing occurs between 
vanes 
— Last fll fraction — fluid completely submerges PMD, different slosh behavior 


Frequencies vs. Fill Fraction 
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Function of pendulum’s length and acceleration 
Monotonic trends 
Frequencies are identical between NTO and MMH 


Frequencies for the two pendulums converge as fill fraction 
increases 


— Sector and annular slosh modes become less distinct as PMD becomes 
fully submerged 


d oN Frequency - Literature Comparison1 


Frequencies vs. Fill Fraction, 


Cassini Paper Frequencies vs. Fill Fraction Comparing to SP-106 Bare Tank 
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+ Left: Cassini paper referenced SP-106 for an analytical equation for 
slosh frequency in a bare tank (cylindrical tank with no PMD) and 
compared it to the frequencies of their two pendulums 

+ Right: Similar trends to Cassini found in Europa pendulum model 
frequencies 

+ Sector and annular slosh modes converge towards bare tank frequency 


as PMD becomes more submerged (fully submerged at 85% fill fraction 
for Europa tank) 
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bs ie 
of Frequency — Literature Comparison2 (ue; 


Frequencies vs. Fill Fraction, 
Comparing to Analytical Sector and Annular Tanks 
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SP-106 references tables (Bauer, 1963) for an analytical equations for sector and 
annular slosh frequency 


Function of acceleration, geometry, and fluid height 
Pendulum frequencies are close (0 analytical equation frequencies 


Differences between analytical and pendulum fits due to: 
— PMD is not exactly a sector/annular tank 


Half-dome bottom approximated as flat bottom — at 25% fill fraction, sloshing fluid is almost 
entirely in the dome 


PMD doesn't include entire height of tank — at 85% fill fraction, PMD is completely submerged 


i, i Damping Ratio 


Damping Ratio vs. Fill Fraction 
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* Monotonic trends 
* Slightly higher damping ratio for higher dynamic 
viscosity (MMH) 


Damping Ratio vs. Fill Fraction 


CFD Fits and Analytical 
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Frequency (rad/s) 
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Mikishev and Dorozhkin found correlation for 
damping in a bare tank 

Function of geometry, acceleration, viscosity, and 
fluid height 

Scales by correction coefficient for domed bottom > 
Pendulum damping within order of magnitude of 
analytical prediction 

Pendulum damping less sensitive to viscosity than 
analytical prediction — viscous vs. drag forces 
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d. " Length and Hinge Location 


Hinge height vs. fill fraction 
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Fill fraction 


+ Origin is top of tank 

+ Pendulum bobs stay within fluid 

+ Monotonic values for pendulum heights 

+ NTO and MMH heights are close but not identical 


NTO 25% fill NTO 50% fill 


Approximate 
enencceees / tank wall 
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PLOTS COMPARING 
PENDULUM MODELS 
AND CFD DATA 


f) NTO 25% Fill Fraction 


Solution Time 2 (s) 


: Volume Fraction of NTO 
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Moment about Y vs Time 
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NTO 50% Fill Fraction NASA 


Solution Time 2 (5) 


CM in X vs Time 
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Force in X vs Time 
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\ NTO 85% Fill Fraction ASA 


Solution Time 2 (s) 
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CM in X vs Time 
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y. » NTO 85% Fill Fraction N 


Force in X vs Time 
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Moment about ¥ vs Time 
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Summary of Parameters NASA 


al NTO (nitrogen tetroxide) MMH (monomethyl hydrazine) 
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Mean Force Error 
from t=0 0.0716 0.075, 0.1055 0.0398 0.0447 0.0679 
Mean Force Error 
from First Peak 0.0241 0.018 0.0775 0.0118 0.0119 0.0518 
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CONCLUSIONS 


Accuracy of Fit NASA 


+ Two-pendulum model can accurately capture either before or 
after first peak 

+ High confidence on frequencies except 85% fill pendulum 2 

+ Moderate confidence on mass, damping, and hinge location 


— Sometimes several sets of parameters could have provided good matching to 
CFD 


— Selected parameters that made physical sense 
+ Model parameters may reflect inaccuracies in CFD 
+ Pendulum model does not scale well for high fluid disturbance 
angles 
+ Damping is actually a function of time and distance traversed by 
moving fluid 
— Pendulum model assumes damping is constant over time 
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+ Small initial fluid displacements: Changes have little 
impact on long-term CFD results 


+ Large initial displacements: behavior differs drastically 


ex Observations to Note NASA 
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* Changing density (NTO vs MMH) only slightly 
changes damping, has little impact on CFD results 


» 
4.4) Areas for Further Investigation N 


* Find literature to support mass fraction parameters 
+ Potentially to capture first peak — add third pendulum 
with damping ratio of one 
+ Validate with more CFD data: 
— Atintermediate fill fractions 


— At different initial fluid offset angles - 5 degree offset is more 
conservative than 15, will be used for deliverable in May 


+ Validate with experiments 


oe 
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